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ABSTRACT 



The analysis of Redshift-Space Distortions (RSD) within galaxy surveys provides con- 
straints on the amplitude of peculiar velocities induced by structure growth, thereby allowing 
tests of General Relativity on extremely large scales. The next generation of galaxy redshift 
surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), and the Euclid ex- 
periment will survey galaxies out to z = 2, over 10, 000-20, 000 deg 2 . In such surveys, 
galaxy pairs with large comoving separation will preferentially have a wide angular sepa- 
ration. In standard plane-parallel theory the displacements of galaxy positions due to RSD are 
assu med to be paralle l fo r all galaxies, but t his assumption will break down for wide-angle 
pairs. ISzapudil d2004l) and Papai & Szapudi (2008) provided a methodology, based on tripo- 
lar spherical harmonics expansion, for computing the redshift-space correlation function for 
all angular galaxy pair separations. In this paper we introduce a new procedure for analysing 
wide-angle effects in numerical simulations. We are able to separate, demonstrate, and fit each 
of the effects described by the wide-angle RSD theory. Our analysis highlights some of the 
nuances of dealing with wide-angle pairs, and shows that the effects are not negligible even 
for relatively small angles. This analysis will help to ensure the full exploitation of future 
surveys for RSD measurements, which are currently confined to pair separations less than 
~ 80 h' 1 Mpc out to z ~ 0.5. 

Key words: cosmological parameters — large scale structure of Universe — cosmology : 
observations — methods : analytical 



1 INTRODUCTION 

Many different mechanisms have been suggested to explain the 
observed late-time acceleration of the e xpansion of the Universe 
dRiess et al.ll998l ; |Perlmutter et aljl999 ) Q Differentiating between 
these options is one of the main challenges facing cosmologists to- 
day. We can try to build up the evidence for different mechanisms 
by examining the evolution of the Universe in two key ways: mea- 
suring the background geometry and measuring structure formation 
within it. 

The geometrical evolution of the Universe can most eas- 
ily be measured using two primary techniques: we can use su- 
pernovae as standard candles or galaxy clustering as a standard 
ruler to make precise measurement of cosmological expansion. Al- 
though supernovae were first to confirm the acce lerated expansion 
of the Universe at hig h statistical significance ( Rie ss et alll 19981 : 
IPerlmutter et alj fl999) following analyses of early galaxy surveys 
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1 Fo r a recent review of dark energy see, e.g., iFrieman. Turner & Hutereil 
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( BfstathiouetalJl990h . supernovae s urveys are no w limited by sys- 
tematic rather than statistical errors ( lKesslerll2009h . The possibility 
of using galaxy clustering to provide a standard ruler has become 
increasingly important since the baryon acoustic peak was detec ted 
dPercival et :afl 1200 ll : ICole et alj 120051 : lEisenstein et~ai1 120051) in 
galaxy power s pectra measured fr om the 2dF Galaxy Redshift Sur- 
vey (2dFG RS:[Colless et alj|2003h and the Sloan Digital Sky Sur- 
vey (SDSS: lYork et alj2OO0l) . Using only the Baryon Acoustic Os- 
cillation (BAO) component of the galaxy clustering signal makes 
constraints robust to non-linear effects, and has already been ex- 
ploi ted to produce intere sti ng constraints on cosmological mod- 
els dPercival etalj|2007a 0; iGaztanaga et alj|2008l : ISanchez et ail 
l2009l;|Percival et alj201fj|T 

Galaxy surveys provide complementary information about 
the build-up of large-scale structure through Redshift Space Dis- 
tortions (RSD). These arise because we do not observe true 
galaxy positions, but instead infer distances from measured red- 
shifts. Coherent comoving galaxy velocities due to the growth 
of structure t herefore lead to measurable anisotropic clustering 
dKaiseJI 19871 ; iHamiltonll 1 998h . RSD have been measured from the 
2dFGRS and SDSS using techniques based on both correlation 
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funct i ons and power-spect ra l Peacock et alj 200 ll; Hawkins et al.l 
20031; iPercivalet all [20041: iPope et al.l 12004: IZehavi et all 1200a 
Okumura et alj2008l : ICabre~ & Gaztanaga 2009), and have recently 
been detected at higher redshift Iguzzo et all 120091 : iBlake et all 
l20ld) . 

In fact, for fundamental reasons, RSD are independent of 
galaxy bias: galaxies act as test particles in the matter flow, so their 
motion is independent of galaxy properties. We can therefore mea- 
sure the matter velocity field at the locations of the galaxies and, 
provided that galaxy positions are representative within the veloc- 
ity field, this gives an unbiased measurement of fog, (mass) where 
/ = d In D/d In a, the logarithmic derivative of the linear growth 
rate, D(z), with respect to the scale factor a, and admass) quan- 
tifies the amplitude of fluctuations in the matter density field. 

Many techniques have been proposed for measu ring RSD, 
includ ing using multipoles of the correlation function ^Hamilton! 
1992), or the power spectrum dPercival & W hite 2009), based on 
the plane-parallel approximation for RSD (IKaiserll 19871) . The ge- 
ometry of the system can be easily incorporated into an analysis 
of the correlation function as the separation of each pair can be 
split into radial and angular components. It is more difficult to per- 
form such a decomposition when measuring the power spectrum, 
although various techniques have been suggested, decomposing the 
dens ity field into a basis of spherical harmonics a nd radial functions 
Ce.g. lFisher et al.ll994l : lHeavens & TavloJl995h . As well as allow- 
ing for the geometry to a survey and RSD within it, we need to 
allow for two effects: 

(i) RSD are degenerate with the angular anisotropy of power 
spectrum caused by the Alcock-Palczinski effect, so measurements 
of growth and geometry from the same survey wil l be correlated, 
and it is sensible to undertake a c ombined analysis dBallinger et al.l 
ll996l : ISimpson & Peacockll2009l). 

(ii) as noted by lHamiltonl fl998). the full RSD operator has extra 
terms compared to the plane-parallel one and results in a redshift- 
space power spectrum which is not diagonal in k. These extra terms 
are important for galaxy pairs separated by wide angle s and need to 
be included in any analysis that uses wide-a ngle data ( Szala v et"al] 
1 19981 : ISzapudill2004 IPapai & Szapudill2008l) . 



In this pap er we concentrate on the second of these effects. Szapudi 
J2004I) and lPapai & Szapudi (2008) proposed a method of comput- 
ing the redshift-space correlation function which does not assume 
a plane-parallel approximation and is applicable to arbitrarily large 
angles. In their treatment the redshift-space correlation function de- 
pends not only on the redshift-space pair separation x and cosine of 
the angle of the pair with respect to the line-of-sight /i, but also on 
the separation angle between two galaxies in a pair 8, which does 
not have to be small. They were able to express the redshift-space 
correlation function explicitly trough a real-space correlation func- 
tion and cosmological parameters. Given the precision that future 
surveys will achieve and the fact that they will cover large fractions 
of the sky, these effects will need to be corrected. 

The wide-angle effects can be subdivided in "purely wide- 
angle" and "mode-coupling" terms: "purely wide-angle" effects 
correct plane-parallel predictions accounting for the fact that the 
separation angle is non-zero, "mode-coupling" terms in addition 
account for the fact that galaxy pairs coherently moved from the 
high-density to low-density regions. These latter terms vanish if the 
initial real-space distribution of galaxies is uniform in distance, and 
both terms vanish in the plane-parallel limit due to the symmetry 
of the system. Both terms are of the same order and mode-coupling 



tends to smooth out features in the power spectrum such as the 
Baryon Acoustic Oscillations (BAO). 

Given the importance and, particularly, the quality of RSD 
constraints expected from forthcoming surv eys such as the Baryo n 
Oscillation Spectroscopic Survey (BOSS: ISchlegel et alj|2009al) . 
surveys resulting from proposals for wide-field mu lti-object spec- 
trograp hs on 4m telescopes such as Big BOSS jSc hlegel et al. 
2009b) and satellite missions such as Euclid dLaureiis et alj|2 009). 
it is timely to revisit this problem, and to test wide-angle RSD the- 
ory. In this paper we use numerical simulations to extensively test 
the dependence of RSD constraints on the angular separation of 
galaxy pairs. 

This paper is organised as follows: in Sec.[2]we briefly review 
the theory of wide-angle RSD; in Sec.[3]we describe a time-efficient 
method of getting a low-noise measurement of wide-angle correla- 
tion function from a mock HV catalog; in Sec. [4] we present our 
results and discuss all the steps necessary to match measured cor- 
relation function with theoretical predictions; in Sec.[5]we conclude 
and discuss how the results will affect real survey measurements. 



2 WIDE-ANGLE REDSHIFT-SPACE DISTORTIONS 

Galaxy positions are measured in a redshift-space, and differ from 
the real-space positions because of the contributions from peculiar 
velocities, 



s(r) 



v r (r)r, 



(1) 



where s is the redshift space position, r is the real space position 
and v r the radial component of peculiar velocity. 

Since the total number of galaxies in real and redshift spaces is 
the same, the number N s (s)d 3 s of galaxies observed in a volume 
element d 3 s of redshift space is related to the real space number 
density N r (r) by: 



N B (s)d 3 s = N r (r)d 3 r. 



(2) 



The observed redshift space galaxy overdensity 8 s at position s can 
be related to the redshift-space selection function N s , 



5 s (s) = 



N s (s) - N a (s) 
N a (s) ' 



(3) 



In contrast, an unbiased estimate of the true galaxy overdensity 
S(r) at position r is given by: 



<T(r) 



N r (r) - N r (r) 
W17) ' 



(4) 



where N r (r) is the expected galaxy distribution in real-space - the 
expected number density of unclustered galaxies at position r given 
the selection criteria of the survey. 

We distinguish here between the measured redshift-space se- 
lection function N 3 {s) and the true selection function N r (r) even 
though the two are the same at linear order. We now assume that the 
expected density is a function of the radial component of the posi- 
tion only, which we denote s in redshift-space and r in real-space. 
Thus the relation between the observed redshift space overdensity 
8" (s) and the true overdensity 8 r (r) is: 

iV s (s)[l + 8 s (s)]s 2 ds = N r (r)[l + 8 r {r)]r 2 dr. (5) 
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Using Eq. {2} we obtain: 

l + S s (s) = [l + <T(r)] ( 1 + 



1+ V - 

r 



dv 
dr 
N(r) 



(6) 



which gives, to linear order: 

5»=<T(r) 

where 

a(r) = 



N(r + vf) 
a(r)v s 



91nr 2 iV s (r) 
dlnr 



(7) 



(8) 



Usually the (1 + term in the Jacobian is omitted because, in the 
linear regime, it gives rise to a 2v/r term, tha t would tend to zero 
at large distances. However, IPapai & Szapudi d2008h have argued 
that, for wide angles, the v/r term is of the same order as the d r v 

term. 

Using the exact Jacobian, and following IPapai & Sza pudi 

(2008), we can express the linear overdensity as: 



d s k 
(2^)3 



so the correlation function is: 



(5 s ( ri )«5 s *(r 2 )) 



rk 



-r-2) 



S(k), (9) 



1+ 3 +Y L2{flk) ~7^ Ll{ '~ nk) 

^ + t + %L 2 {r 2 k)+ l -^-L 1 {r 2 k) 
6 6 r 2 fc 



(10) 



where Li are the Legendre polynomials of order i and P(k) is the 
real-space matter power spectrum. The third terms in the brackets 
are the ones responsible for the wide-angle effects, while the fourth 
terms are the ones responsible for the mode-coupling. Note that the 
?-i and r 2 terms in the denominator depend on the angular separa- 
tion of the galaxies. Setting <f>i as the angle between the vector to 
the first galaxy in a pair t%, x to be the vector connecting galax- 
ies in a pair, and (f>2 to be the angle between vector to the second 
galaxy in a pair r 2 and x, we can use the sine rule (see Fig.[T} to 
express n and r 2 as 



71 



sin(< 



T2 = 



(11) 



(12) 



Tripolar spherical harmonics are the most natural basis for 
the expansion of a fun ction that depends on three direc tions 
l Varshalovich et alj 19881) so, as suggested in Szapudi (2004) and 
IPapai & Szapudild2008l) . we expand Eq. dlOt using a subset of them 
that are proportional to the zero angular momentum: 



Se 1 e 2 l(ri,f'2,x) 



(13) 



£i 

im 



in 2 



77! 



imi {fl)Ce 2 m 2 (r2)Cirn{x) , 



where Ce m are normalized spherical functions, multiplied by the 
3j Wigner symbol. The redshift space correlation function is then 
written as: 

£ s (fi, f 2 ,x)= ^2 B lll2 \x,(j> 1 ,<j>2)St 1 i 2 t(r 1 ,r2,x), (14) 
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Figure 1. The coordinate system adopted for the Mangle formed by the 
observer O, and galaxies Gi and G2- 



where B lll2l (x, cf>i, 4>2) are a series of coefficients that depend on 
/, gi(<pi) and £j ( x), and they can be di vided into wide -angle com- 
ponents (given in lSzalav et alJfl99^ & ISzap udi 2004) that do not 
depend o n the third term in Eq. l[9j, and mode-coupling components 
(given in IPapai & Szapud i 200 8); the plane-parallel approximation 
emerges as a limit wh en f i = V2- 

As described in lHamilton & Culhand ( fl99r3) . for any pair of 
galaxies we can fully analyse the problem within a plane formed 
by the two galaxies and the observer; the effect of redshift-space 
distortions only depends on the geometry within this plane. In the 
coordinate system described above, the vectors connecting the ob- 
server with the galaxies are ri and r?, and the separation vector 
between the galaxies is x; 4>i are the angles between ri and x. In 
the plane parallel approximation cj>i = cj>2- This description of the 
triangle is shown in Fig.Q] 

At this point we still have the freedom to choose a line-of-sight 
to the galaxy pair, we follow the standard approach and choose the 
line of sight as a direction bisecting the angle formed by the two 
galaxies, which we assume is given by 29. The angle between the 
angle bisector and x is denoted 0, and we define [i = cos(0). We 
therefore have that 4> = \{4>i + 02 ), and 9 = ^(4>2 — (j>i). This is 
shown in Fig.Q] 

With this choice of coordinate system the redshift-space cor- 
relation function can be expanded as: 



e{x,o,<t>) = J2 Cab ^ L ^ cos ^ L ^ cos ^' (15) 

ab 

= ^ a ni n 2 cos(ni^i) cos(n 2 <fe) 

+b„ 1Tl2 sin(ni0i)sin(n 2 <?!>2), (16) 



with coefficients given in IPapai & Szapudi fc008h . Eq. (B) reduces 
to the plane-parallel result ( Ha miltonlll992h in the limit of 9 — > 0. 
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3 MEASURING WIDE-ANGLE CORRELATION 
FUNCTIONS 

A standard method for analysing a N-body simulation to test wide- 
angle redshift-space distortions would be to: 

(i) locate an observer within the output, 

(ii) translate all galaxies from real into redshift-space based on 
this observer, 

(iii) sample from these galaxies based on desired radial distribu- 
tion, 

(iv) split pairs into bins in <f>, 9, x and counts pairs, 

(v) estimate the correlation function. 

This approach mimics that of an actual survey analysis, creating a 
mock galaxy catalogue. However, if we want to measure the effect 
of redshift-space distortions for galaxy pairs with a particular angu- 
lar separation, the method is not optimal as only a small fraction of 
the galaxy pairs analysed will have this angular separation. In addi- 
tion, we have to perform the full procedure for every radial galaxy 
distribution that we wish to analyse, and this distribution will limit 
the pair-density that can be selected from the simulation. In order 
to rapidly increase the signal-to-noise, we adopt a different proce- 
dure, allowing the origin to move so that each galaxy pair can be 
analysed as if it was observed with the required angular separation. 
This procedure can be summarised as follows: 

(i) decide on the value of 9 for which we wish to analyse pairs, 

(ii) take each galaxy pair from the simulation with real-space 
separation < i? ma x h -1 Mpc, 

(iii) for each pair randomly choose fi £ [— cos(6 l ), cos(#)], 

(iv) choose the location of the origin giving this /i and 9, 

(v) move galaxies according to their expected redshift-space 
distortion, 

(vi) weight the pair by a function of /i, x to match desired dis- 
tribution, 

(vii) split pairs into bins in <j), x and counts pairs, 

(viii) estimate the correlation function. 

The added complexity of including redshift-space distortions on a 
pair-by-pair basis is outweighed by the ability to obtain more pairs 
at the desired angular separation. Note that /i in step (iii) is con- 
strained and cannot have any value within [—1,1] because \<j}\ > \9\ 
results in geometrically impossible triangles (see, Fig.[TJ. The dis- 
tributions of triangles in /i and x for a galaxy distribution with a 
power law selection in r can be calculated analytically, as we now 
demonstrate. 



where dNp lL1T p(x, /j.) is the number of pairs with separation x ± dx 
and cosine of the angle with the line of sight in the /i ± d/i interval. 
If we pick galaxy pairs separated by a fixed opening angle from the 
uniform spatial distribution of galaxies (step (i) in the procedure 
described above), for a large enough angle the distribution of pairs 
will not follow Eq. For galaxy pairs with a fixed half-angle 
9 (see, Fig. QJ, the likelihood of finding one galaxy at position ri 
and another at v% is P(ri) oc r\dri, P{r2) oc r 2 dr2- Since the 
two likelihoods are independent, the joint probability of finding that 
pair of galaxies is: 



dNg ml (r\,r2) oc r 1 r 2 dr\dr2 



(18) 



Using Eqs. i ll It and i l 1 2b we can rewrite Eq. J 1 8b in terms of cf> and 

x as: 



dNg alI (x, fi) oc x sin(</> + i 



' — 9) dxd(f>. 



(19) 



After completing steps (ii) and (iii) in the previous subsection, we 
will have a distribution of galaxy pairs that is uniform in dfi — 
s'm((f>)d(f> and scales as x 2 dx with separation. This distribution of 
galaxy pairs does not correspond to a = 2, If we want to compare 
our data with the correlation function computed from Eq. J 1 5 1 > for 
q = 2, we have to weight our galaxy pairs (step (iv) in the previous 
subsection) by an additional factor of x 3 sm((f> — 9) 2 sin(<^ + 9) 2 
to get a distribution given by Eq. H9\ . 

This procedure can be applied to the case of arbitrary a. If, for 
example, the galaxy number density scales as a power-law n(r) = 
r~ , then a — 2 — TV, and the distribution of galaxy pairs with a 
fixed opening half-angle drawn from this distribution will be: 



dNS 



r (x,fi) oc x 2a+1 sin( 



< + 9) a sm{<t>-9) a dxd<t>. (20) 



We therefore weight all pairs recovered from the simulation with a 
weight given by: 



W(x,fi) = 



+ 9) a sxa(<pr - 9) c 



(21) 



where x s is the separation of the galaxies in redshift-space, (f> r is 
the real-space value of <j> for the chosen galaxy pair and <f> s is the 
redshift-space value. We divide the real-space distribution by the 
redshift-space equivalent in order to normalise the weights to give 
no net change in the expected pair distribution. 



3.1 Real-space distribution of fi, x 

The redshift-space correlation function in Eq. d!5t depends on the 
selection function a(r) = 2 + d\n[n(r)]/dr, where n(r) is the 
number density of galaxies in real space and r is the position of the 
galaxy from the observer as in Section[2] If the galaxy distribution 
is uniform (a(r) = 2, n(r) = n, both independent of r), then the 
probability of finding a galaxy in some region is proportional to the 
volume of that region. Therefore if we randomly pick a galaxy in 
the catalog, there will be on average nx 2 dxdfi galaxies in a small 
volume element dxdfi which is distance x ± dx away from that 
galaxy within an angular slice /i ± d[i. In the plane-parallel ap- 
proximation, where all galaxies are assumed to lie along the same 
direction from the observer, the distribution of galaxy pairs will 
scale as 

dN^ iT (x, n) oc x 2 dxdfi, (17) 



3.2 Redshift-space distribution of fi, x 

When transformed into redshift space, the real-space galaxy distri- 
bution is "washed out" by the random component of galaxy veloci- 
ties, so we need to convolve with the random velocity distribution. 
This would normally not be a problem as we would estimate the 
galaxy distribution directly in redshift-space. However, using the 
above procedure we have set the real-space distribution of galax- 
ies, so we need to take care when modelling the redshift-space dis- 
tribution. To correct for this effect, we first estimate the velocity 
dispersion a v from the catalog, assuming that the random veloci- 
ties are drawn from a Gaussian distribution. We then convolve the 
initial distribution in real space with this Gaussian. For example, 
for n(r) oc r~ N , the probability of finding a galaxy at distance r 
away from the observer is P(r) oc r 2 ~ N — r a . In the unclustered 
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5 (x) for B = 0.355 
= 2 data without gaussian convolution 
= 2 data with gaussian convolution 
= 2 full mode coupling theory 



Figure 2. The monopole correlation function calculated from the HV mock 
catalogues designed to mimic a galaxy distribution with 9 = 0.355, and 
a = 2. Open circles (with Icr errors) were calculated using the redshift- 
space model for the galaxy distribution. Open diamonds used the real-space 
model. The solid line shows the prediction from Section|2] 



redshift space this will transform into: 



P(s) 



n /2tFctJ 



(ct„) 



a + 1 



2 



v/2 
a+1 



a + 2 



t-Fi 



_a 3 £_\ 

' 2 ' 2a% J 



iFi 



-a 1 s 2 \ 
~ , 2 , ~2^J 

where F is a hypergeometric function. 

For a specific case of a = 2, Eq. i22l results in: 



(22) 



P( S )oc(cr^+ S 2 )(l + Erf (-g- 



2 

— a v se 

7T 



3 



(23) 

where Erf(r) is an error function. The distribution of galaxy pairs 
for a = 2 will be: 



A£S oc P{ Sl )P{sz)x s dx s d<t>, 



(24) 



where P(sj) is given by Eq. J23t and si and S2 can be expressed 
in terms of s and <j> via Eqs. (TTJ and d 1 2b . after the substitution of 
x with x s . 

When modelling the theoretical correlation function, 
lHamilton & Culhand dl996l) argue that the difference between 
n(s) and n(r) is small. While this is true when creating a model 
correlation function, it is not true, and we need to use the n(s) 
instead of n(r), when we model the data. This is shown in Fig. [2] 
where we plot £o(x) calculated for a — 2 from a mock catalog, 
assuming that the galaxy distribution follows the one expected 
from either the real-space or redshift-space calculation. See 
Section[4]for more details of the calculation. 



3.3 Pairs close to the origin 

Galaxies that are close to the origin can cause problems as, for such 
galaxies, the redshift-space displacement can be larger than the dis- 
tance to the galaxy. This problem is exacerbated because we do not 
include the velocity of the observer in our calculation. In extreme 
situations, naively applying the expected redshift-space distortion 
would place the galaxy on the opposite side of an observer. In or- 
der to avoid such problems, we only include galaxy pairs where 



G2 \ 




( G1 


^JC>^1^\ 0=0.355 y^J^CI^ 




0=0.175 




Ol 




Ol 




0=0.11875 





Figure 3. Schematic representation showing the positions at which the ori- 
gin could be placed for different values of 9. The shaded circles around 
each galaxy are exclusion zones within which the redshift-space distortions 
is comparable with the distance to the galaxy. Ol and 02 mark two pos- 
sible positions for an origin that would give angular galaxy separation of 
29 = 0.71. 



both galaxies are more than 5a v away from the origin, where a v 
is the ID velocity dispersion of the galaxy population. For the HV 
simulation cr„, = 3.9 h~ Mpc. Fig.[5]shows, for galaxies Gi 
and G*2, the loci of positions at which the origin could be placed for 
fixed 8. The circles mark the exclusion zones. If one of the galaxies 
in a pair is inside the exclusion zone in real space, we do not include 
that pair when we estimate the correlation function. This exclusion 
is tracked when we calculate the expected galaxy distribution. 



3.4 The Hubble Volume Simulation 

We apply the procedure outlined in Sec. [3] to analyse wide-angle 
redshift-sp ace distortions with in the ACDM Hubble Volume (HV) 
simulation fevrard et al]|2002l) . The ACDM HV simulation, cover- 
ing a (3000 h" 1 Mpc) 3 box, assumes a cosmological model with 

n m = 0.3, ficDM = 0.25, n b = 0.05, n A = 0.7, h = 70, 

erg = 0.9, & n s — 1. We do not apply any galaxy bias, simply 
Poisson sampling the matter particles to give our "galaxy distribu- 
tion"; the inclusion of a bias model would not alter the conclusions 
of this work. We use the periodic nature of the numerical simu- 
lation to eliminate boundaries from our pair counts. This means 
that, by using the above weighting scheme and allowing for the re- 
moval of galaxies close to the origin, the expected number of galaxy 
pairs in the absence of clustering RR can be calculated analyti- 
cally. We can therefore use the natural estimator 1 + £ = DD/RR 
dLandv & Szalav|[l993h , where DD is the measured number of 
galaxy-galaxy pairs. 



4 RESULTS 

We have performed 100 runs, each based on a sample of 10 6 galax- 
ies drawn from the z — output of the HV ACDM simulation. For 
each unique pair of galaxies within the sample we have selected two 
locations for the origin, one where the galaxies su bject a separation 
angle of 9 — 0.355, chosen to match figure 1 of Papai & Szapudi 
(2008), and one origin where 6 — 0.71, twice this angular sepa- 
ration, to emphasise the wide angle effects on pairs separated by a 
very wide angle. 

We include all pairs with redshift-space separation x s < 
200 h~ x Mpc, and each pair was weighted as described in Sec- 
tion !3.2l based on its real-space separation and angle to the line of 
sight. Given that, as one can see from Eq. dlOt , the mode-coupling 
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o data 

wide angles + mode coupling theory 

— pure wide angles theory 
plane parallel theory 



Figure 4. Correlation function as a function of scale for fixed 6 = 0.355 
and for a bin c e ntered at = 1.03, designed to match figure 1 of 
Papai & Szapudi 12008). The plane-parallel model is shown by the dotted 
line, the pure wide-angles model by the dashed line and the full mode- 
coupling prediction with a = 2 by the solid line. We expect the match 
between data and models to be good at scales larger than 40 ft -1 Mpc, 
because we didn't model the small scale non linearities. 



terms depend strongly on the radial galaxy distribution, for both 
the aperture angles we selected samples with different radial galaxy 
distributions, in order to test our methodology in different cases; we 
chose a = 0, 0.5, 2, 4: 

• a = corresponds to the case of a galaxy distribution that has 
equal density in radial bins of equal width. In this case there are no 
mode-coupling terms in Eq. dlO) . so this corresponds to the "pure 
wide angles" case. 

• a = 0.5 corresponds to a final real-space distribution of 
galaxy pair separations with probability density function propor- 
tional to x 2 dx. This is the same distribution obtained by randomly 
sampling pairs from the simulation, so the x factor drops out of 
Eq. \2\\ . Note that, for clarity, we do not present results from this 
value of a as they are very similar to, and overlap results for a = 0. 
The match between data and theory has the same quality as for the 
other values of a. 

• a = 2 matches a distribution of galaxies that are uniformly 
distributed in volume, so equal volumes contain equal numbers of 
galaxies. Most planned surveys aim to observe galaxies with this 
radial distribution. Note that, for this galaxy distribution, the pair 
distribution goes as x 5 dx. 

• a = 4 represents a steeper radial distribution where more 
galaxies are found at larger distances. This will increase the ef- 
fect of the mode-coupling terms in the wide-angle redshift-space 
distortion formulae (Eq.l21t. 

We have calculated errors by comparing outputs from 100 sub- 
samples, that cover the same HV volume and so only contain the 
shot noise element. These errors will therefore underestimate the 
true error, because they do not fully include the cosmic variance 
component. However, the HV volume is (3000 ft -1 Mpc) 3 , and we 
only consider pairs with x < 200 h~ x Mpc, so we expect the shot 
noise to dominate the error budget. Even so, it is worth pointing out 
that our primary aim is to consider deviations from plane-parallel 
theory and that our proposed methodology works to match data 
with the full mode-coupling theory: the size of the errors is unim- 
portant, provided they are far smaller than the differences between 
theories. 

Fig-0shows the correlation function, calculated within a nar- 



8 (rad) 

Figure 5. Ratio of theoretical full mode-coupling correlation function to 
plane-parallel correlation function, computed for a = 2 at fixed value of 
r = 79 h —1 Mpc; the points mark the values of 8 used in our simulations. 



row bin in n for 6 = 0.3 5 5 and a = 2; this is equivalent to fig- 
ure 1 of IPapai & Szapudil d2008h . We are able to fit the theoreti- 
cal correlation function to the estimate from HV simulations for 
scales larger then 40 ft" 1 Mpc (we do not expect to match per- 
fectly data with theory on smaller scales because we didn't model 
non-linearities). Looking at plane-parallel, pure wide-angle and full 
mode-coupling theories, it is clear that only the full mode-coupling 
theory provides a good fit to the data. 

Fig. [5] shows the ratio between the correlation function pre- 
dicted by the full mode-coupling theory and that for the plane- 
parallel case, as a function of 6, for a fixed value of r = 
79 h~ Mpc, with a = 2; from this plot we can see that even for 
small angles there is a non negligible difference. 

In order to analyse RSD, the measured correlation function is 
usually decompo sed into Legend re momenta, which contain all of 
the RSD signal (Hamilton 1998). Some combination of these mo- 
menta is then used to constrain cosmological parameters through 
their effects on the growth of structure. To see how wide angle 
effects would modify these measurements we estimate first three 
even Legendre momenta of the correlation function from HV mock 
catalog. We measure: 



Y,DD(x,^)L e (^) 



(25) 



for I = 0, 2, 4. Adopting the plane-parallel philosophy ^Hamilton! 
1 19921) . one might be tempted to interpret these functions as: 



£o PP (x) 

£ PP M 



1 + -0 + -0* ) £(x), 



35 



3 + -p)Z(x), 



(26) 
(27) 
(28) 



where f (a;) is a real space correlation function and ft = f/b. In 
fact the functions in Eq. l |25t are given by: 



Wr(ji)?(x,0,n)Lt(jt)dti, 



(29) 



where W T {fj) is a weight function that appears because of geo- 
metrical constraints and because the distribution of fi will not be 
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uniform. We do not correct for this weight in our measurement 
from simulations, and apply it to the theoretical calculations that 
we compare against. 

Fig.[6]presents the main results of this work, showing the ap- 
plication of our methodology to sample and analyse data, and com- 
paring against the full mode-coupling predictions. We plot £o(x), 
&{x), for both the aperture angles and for radial galaxy den- 
sity distributions corresponding to a = 0, a = 2 and a = 4; as 
we can see, the contribution of the mode coupling terms increase 
with a, and as the radial galaxy distribution steepens, more galax- 
ies are moved nearer by the RSD, leading to an increase in the 
small-scale correlation function. Momenta of the correlation func- 
tion computed with the full mode coupling theory match remark- 
ably well the data analysed as explained in Sec. [3] We also plot 
the plane-parallel prediction as a comparison, and this demonstrate 
how badly this approximation fails for galaxy pairs with wide angu- 
lar separation. The only parts where our methodology do not match 
the mode coupling theory are where the correlation multipoles are 
very small (around 10 -4 ) or for small scales, where non-linearities 
become important; we also have to take in account that, for the 
steeper galaxy distributions, we are upweighting results from the 
same number of pairs as the more shallow galaxy distributions so, 
although the signal is stronger, the relative error will be the same. 



5 DISCUSSION AND CONCLUSIONS 

Redshift-Space Distortion Analyses are a powerful tool for cosmol- 
ogy, but incoming data need to be analysed very carefully in order 
to fully extract all available information. Until now RSD analyses 
have concentrated on using the plane-parallel approximation (with 
a few exceptions including lOkumura et al. 2008; Pop e et al ] |2004 
Matsubar l2004h ; this is almost correct if the survey is narrow, but 
when galaxy surveys with wide field of view data will be available, 
there will be a consistent number of galaxy pairs separated by wide 
angles and in this case the plane-parallel approximation fails. 

This approximation arises from not taking in account some of 
the terms in the Jacobian relating redshift- to real-space; in this pa- 
per we considered its exact expression, and this causes additional, 
non-diagonal terms in the correlation function. We then use the ex- 
pansion of the correlation f unction in a base of tri polar spherical 
harmonics, as suggeste d in Papai & S zapudi 20 0j|; follo wing the 
formalism developed in lSzapudill2004l ; |Papai & S zapudi 200j|, we 
make predictions for the momenta of the correlation function for 
galaxy pairs at fixed angular separation, and we test them against 
data from the Hubble Volume Simulation. 

In order to do this we have had to introduce a new method- 
ology: rather than creating a single sample of galaxies in redshift- 
space from which we can count pairs, we have instead dynamically 
applied RSD on a pair-by-pair basis, choosing an origin for each. 
By including weighting functions in x and fi we can match results 
from the more traditional approach. This allows us to only consider 
galaxy pairs at particular values of 9. 

To show both the correctness of our methodology and the de- 
viation from the plane-parallel situation we tested galaxies sepa- 
rated by two different fixed values of 9; the mode coupling terms 
are also strongly dependent on the galaxy radial distribution, so we 
have tested simulations with 4 different number density distribu- 
tions, for both values of 0. 

We show that taking in account the wide angl e and mode cou- 
pling terms (that are of the same order, as stated in lPapai & Szapudil 
2008) give a clear deviation from the plane-parallel theory; using 



the exact theory and our methodology, we can match the results of 
simulations, the agreement between data and theory being remark- 
able, especially considering how crude our modelling of the HV 
simulations is (we use a measurement of the 3D real-space cor- 
relation function as our baseline model, and we do not include a 
correction for fingers-of-god type effects). For RSD measurements 
made within radial bins, it will be vital to match the theory to the 
exact distribution of galaxies observed. 

In a measurement of RSD from real data, the final result will 
be a weighted average over different opening angles, to take ac- 
count for the fraction of galaxy pairs separated by a 9 angle, so in 
wide surveys one has to discard the plane-parallel approximation, 
or face losing a considerable amount of information. For this reason 
our methodology will be particularly useful for incoming and fu- 
ture redshift Surv eys, such as the Baryo n Oscillation S pectroscopic 



Survey (BOSS; ISchlegel et al] l2009ah. BigBOSS dSchlegel et ail 
2009b), and Euclid dLaureiis et alj |2009); the importance of this 



methodology for curre nt and future experiments will b e considered 
in a subsequent paper dSamushia et aljin prepara tion). 



6 ACKNOWLEDGEMENTS 

AR is grateful for the support from a UK Science and Technol- 
ogy Facilities Research Council (STFC) PhD studentship. WJP is 
grateful for support from the European Research Council and from 
the Leverhulme Trust and STFC. LS acknowledges support from 
European Research Council, Georgian National Science Founda- 
tion grant ST08/4-442 and SNSF (SCOPES grant No 128040). AR 
would like to thank G. W. Pettinari for the useful discussions. Sim- 
ulated data was calculated and analysed using the COSMOS Altix 
3700 supercomputer, a UK-CCC facility supported by HEFCE and 
STFC in cooperation with CGI/Intel. 



Peacock J. E., Heavens A. E, 1996, MNRAS, 



REFERENCES 

Ballinger W. E., 
282, 877 

Blake C, et al., 2010, MNRAS accepted, larXiv:1003.57211 

Cabre A., Gaztanaga E., 2009, MNRAS, 396, 1119 

Cole S., et al., 2005, MNRAS, 362, 505 

Colless M., et al., 2003, |astro-ph/0306581| 

Eisenstein D.J., et al., 2005, ApJ, 633, 560 

Evrard A.E., et al., 2002 ApJ, 573, 7 

Efstathiou G, Sutherland W.J., Maddox S.J., 1990, Nature, 348, 
705 

Fisher K.B., Scharf C.A., Lahav O, 1994, MNRAS, 266, 219 

Frieman J., Turner M., Huterer D., 2008, Ann. Rev. Astron. As- 
trophys., 46, 385 

Gaztanaga E., Cabre A., Hui L., 2008, | arXiv:0807.355T1 

Guzzo, L., 2009, Nature, 541, 2008 

Hamilton A.J.S., 1992, ApJ, 385, L5 

Hamilton A.J.S., Culhane M., 1996, MNRAS 278, 73 

Hamilton A.J.S., "Linear redshift distortions: A review", in "The 
Evolving Universe", ed. D. Hamilton, pp. 185-275 (Kluwer Aca- 
demic, 1998) |astro-ph/9708T02) 

Hawkins E. et al., 2003, MNRAS, 346, 78 

Heavens A.F., Taylor A.N., 1995, MNRAS, 275, 483 

Kaiser N., 1987, MNRAS, 227, 1 

Kessler R., et al., 2009, ApJS, 185, 32 

Landy S. D., Szalay A. S., 1993, ApJ, 412, 64 



© 2010 RAS, MNRAS 000,[T]l9] 



8 A. Raccanelli et ah 



a (x) for = 0.355 
a = data © a = 4 data 

a = theory a = 4 theory 

a = 2 data plane-parallel theory 

a - 2 theory 




\m v 



20 30 40 50 60 70 SO 90 100 110 120 13(1 140 150 



yx) for 9 = 0.355 
a = data o a = 4 data 

a - theory a - 4 theory 

a = 2 data plane-parallel theory 

a = 2 theory 




30 40 50 60 70 SO 90 100 110 120 130 140 




^(x) for = 0.355 
a a = data o a = 4 data 

a - theory a - 4 theory 

o a = 2 data plane-parallel theory 

— a = 2 theory 



^(x) for - 0.71 
a = data © a = 4 data 

a = theory a = 4 theory 

a = 2 data — plane parallel theory 
■ a = 2 theory 



20 30 40 50 60 70 SO 90 100 110 120 130 140 150 

X 

^(x) for 0-0.71 
A a = data © a = 4 data 

a = theory a = 4 theory 

plane parallel theory 




20 30 40 .SO 60 70 SO 90 100 110 120 130 140 150 



£ 4 (x) for 9 = 0.71 



a = 4 data 

a = 4 theory 

plane parallel theory 



20 30 40 50 60 70 



90 100 110 120 130 140 150 




20 30 40 50 60 70 



00 100 110 120 130 140 150 



Figure 6. !;q(x), £2(2')" ^i( x ) f° r angular galaxy pair separations of 9 = 0.355 (left column) and 8 = 0.71 (right column) for different power law radial 
galaxy density distributions; lines are the expected values from the theory of Papai & Szapudi 12008) with wide angle and mode-coupling effects, with a = 
(dotted), a = 2 (solid), and a = 4 (dashed); dot-dashed lines are plane parallel predictions. Symbols were measured from the HV mocks, with a = 
(triangles), a = 2 (circles), and a = 4 (diamonds). 



Laureijs R., et al., 2009, "Euclid Assessment Study Report for the 

ESA Cosmic Visions" 
Matsubara T., 2004, ApJ, 615, 573 

Okumura T., Matsubara T., Eisenstein D.J., Kayo I., Hikage C, 

Szalay A.S., Schneider D.P., 2008, ApJ, 676, 889 
Papai P., Szapudi I., 2008, MNRAS, 389, 292 
Peacock J.A., et al., 2001, Nature, 410, 169 
Percival W.J., et al., 2001, MNRAS, 327, 1297 
Percival W.J., et al., 2004, MNRAS, 353, 1201 
Percival W.J., et al., 2007a, ApJ, 657, 645 

Percival W.J., Cole S., Eisenstein D., Nichol R., Peacock J.A., 

Pope A., Szalay A., 2007b, MNRAS, 381, 1053 
Percival W.J., White M„ 2009, MNRAS, 393, 297 



Percival W.J., et al., 2010, MNRAS, 401, 2148 

Perlmutter S., et al., 1999, ApJ, 517, 565 

Pope A.C., et al., 2004, ApJ, 607, 655 

Riess A.G., et al., 1998, AJ, 116, 1009 

Samushia L., Raccanelli A., Percival W.J., in preparation 

Sanchez A.G., Crocce M., Cabre A., Baugh CM., Gaztanaga E., 

2009, MNRAS, 400, 1643 
Schlegel D., White M., Eisenstein D.J., 2009a, larXiv:0902 .4680 1 
Schlegel D., et al., 2009b, larXiv:0 904.0468 1 
Simpson, E, Peacock, J. A., 2009, larXiv:0910.38341 
Szalay A.S., Matsubara T., Landy S.D., 1998, ApJ, 498, LI 
Szapudi I., 2004, PhRvD, 70, 083536 

Varshalovich D. A., Moskalev A. N„ & Khershonski V. K., 1988, 



© 2010 RAS, MNRAS 000,[T]l2] 



Redshift-Space Distortions at Wide Angles 

Quantum Theory of Angular Momentum (Singapore: World Sci- 
entific) 

York, D.G., et al., 2000, AJ, 120, 1579 
Zehavi, I., et al., 2005, ApJ, 630, 1 



© 2010 RAS, MNRAS 000,[T]l9] 



